function [coeff,se] = cov_hac(y,x)

% coeff is cov(y,x)
% se is the standard error
% x needs to be one vector
[coeff_temp,se_temp] = beta_hac(y,x);
coeff = coeff_temp.*var(x);
se = se_temp.*var(x);

end

